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We show that current locahzation in polycrystalline varistors occurs on paths which are, usually, 
in the universality class of the directed polymer in a random medium. We also show that in ceramic 
superconductors, voltage localizes on a surface which maps to an Ising domain wall. The emergence 
of these manifolds is explained and their structure is illustrated using direct solution of non-linear 
resistor networks. 



Low energy paths and surfaces in random systems are 
of broad importance in statistical mechanics and in ma- 
terials' theory. Examples include the morphology of flux 
lines in the presence of random pinning j^] , magnetic do- 
main walls in random-bond Ising magnets [^J^l, dielec- 
tric breakdown paths and fracture surfaces Q . The 
directed polymer in a random medium has received an 
enormous amount of interest, both due to its intrinsic 
importance and due to its relation with growth processes 
. Domain walls in random-bond Ising models have also 
received a good deal of attention, in fact the directed 
polymer in a random medium was first used to model 
domain walls in two-dimensional random-bond magnets 
1^ . More recently it has been realised that random path 
and surface problems map to classic problems in com- 
puter science In particular, the shortest-path prob- 
lem is related to the directed polymer in a random me- 
dia ]lO| and the minimum- cut problem is related to do- 
main walls in random Ising magnets [ pl| . The shortest 
path and minimum cut problems have polynomial time 
algorithms, and this has enabled physicists to solve path 
roblems with overhangs (i.e. not just the directed case) 
and domain wall problems in three dimensions [pl , 
both of which were previously considered difficult. Many 
more connections between statistical physics and solv- 
able combinatorial optimization problems have recently 
proven profitable Here we show that the shortest 
path and minimum-cut problems also emerge in the study 
of the onset of current or voltage in non-linear resistor 
networks and in their applications to varistors and su- 
perconductors. This enables a detailed characterization 
of the manifolds on which current or voltage localize near 
critical thresholds. 

Non-linear flow in random networks is germane to 
problems ranging from non-Newtonian fluid flow in 
porous media u^, to the onset of flow in varistors 
p^-p7[, and to the onset of voltage in grain-boundary 
limited superconductors pq-|20|. The emergence of spe- 



cial flow paths in the fluid and varistor cases is well- 
established and network models with realistic non-linear 
flow laws have been constructed [p^-|l5| . Recently lattice 
models have also been used to model the behavior of poly- 
crystalline superconductors |20|. The networks which 
model these materials have highly non-linear current- 
voltage (or pressure-flow) behavior on each bond of a 
graph. In the case of diodes or varistors each bond has 
a critical voltage, wj, which varies from bond to bond 



jl| , p2 141^]. The onset of current flow occurs on spe- 
cial paths and the identification and characterisation of 
these paths is a central issue in the analysis of varistors. 
In the superconductor case, each bond has a local critical 
current, i^, below which the bond carries current but has 
zero voltage. Above this threshold current, a local volt- 
age must be applied to the bond in order to increase the 
current further. The onset of voltage occurs on a surface 
pof which we show maps to a minimum cut through the 
network. 

Firstly we define the shortest path and minimum-cut 
problems [^,^. These problems (and the algorithms 
which solve them) are very basic to algorithmic systems 
in computer science. They are defined on a graph where 
each edge is assigned a cost (shortest-path problem) or a 
capacity (minimum-cut problem) . Consider then a graph 
composed of nodes and edges, where each edge has a cost, 
Cfe. The shortest path between any two sites, s and t, in 
the graph is simply the path on which the sum of the 
edge costs. 



Cst = minp(y^Cfc), 



(1) 



fecP 



on the path is smallest. In physics language the shortest 
path is the lowest energy path. It is clear that this is the 
same as the problem of a directed polymer in a random 
media, provided the path is forced to be directed. How- 
ever if overhangs are allowed this problem can also be 
solved efficiently using, for example, Dijkstra's method. 
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It turns out that overhangs do not change the universal- 
ity class of the problem unless disorder is strong 24|. 

The minimum-cut problem is defined on a graph 
consisting of nodes and edges in which each edge has a 
flow capacity, Uk- This is the maximum amount of flow 
that the edge can carry. There is no energy cost in the 
computer science definition of this problem, there is just 
the limiting capacity, Uk (However in the end some of 
the physics applications associate the capacity with an 
energy). Now consider the amount of flow that can be 
carried between two sites, s and t. A key theorem states 
that the maximum flow between s and t is equivalent to 
the capacity of the minimum cut between s and t . 
The minimum cut is the lowest capacity surface of sepa- 
ration in which s is on one side of the cut and t is on the 
other. The capacity of the minimum cut, C/st, is the sum 
of the capacities on this surface of separation, S*, 



Ust = ming(y^ 



Uk) 



(2) 



kcS 



If we associate the bond capacities with the exchange 
constants in a random-bond Ising magnet, then the min- 
imum cut maps to the lowest energy domain wall. Some- 
what surprisingly, the minimum cut can be found effi- 
ciently (in almost linear time in the number of edges) 
and this has enabled detailed study of domain walls in 
three dimensions pdlJsot] . 

Now we consider in more detail the emergence of 
the shortest path and the minimum cut in non-linear 
random resistor networks. To be concrete, we con- 
sider square networks where each bond has a monotonic 
(single-valued) non-linear V-I behavior. The case of cur- 
rent onset at a critical voltage is applicable to varistors 
[^(Fig. 1) and diodes and also to Bingham plastics 
(replace voltage by pressure and current flow by fluid 
flow so that flow initiates at a critical pressure). The 
case of a voltage onset at a critical current is applicable 
to ideal josephson junctions (Fig. 2) and also to ideal 
flux-flow laws pO| . Ideal varistor behavior is mathemati- 
cally described by, i{v) — [io + ^^^^)9{v — Vc), while ideal 
Josephson junction or flux flow behavior is described by, 
v[i) — (wo + ^^)^(* — ic), where 9{x) is the step func- 
tion which is zero for a; < and one for a; > 0. Finding 
the current or fluid flow through these non-linear net- 
works seems difficult due to the non-smooth behavior 
at the critical current or voltage. However a variety of 
regularization methods are available to handle these dif- 
ficulties. The simplest is to replace the step function 
by, e{x) = a;/(^2 _,_ 2,2)1/2 ^g^jjg ^jjg Yimit ^ ^ 0. 

This form has nice analytic properties which assist in the 
numerical analysis, so we have used this regularization 
in our codes. A key observation is that these problems 
have a related cost function defined by, cost = v{i')di' , 
which is convex provided the V-I law on each bond is 
monotonically increasing. With suitable regularization 
this is true and solving these non-linear resistor networks 
then reduces to a convex optimization problem with the 



constraint of flow conservation at each node. We have 
used this framework to develop efficient codes for this 
class of problem. Here we discuss applications to varis- 
tors and superconductors. 

In the varistor case (see Fig. 1) each bond has a ran- 
dom onset voltage, v^, and the whole network has the 
onset voltage, Vc- At 14, current only flows on the short- 
est path through the network, as we confirmed numeri- 
cally by finding the path on which Vc = J2keP ''^c is min- 
imal (using Dijkstra's method). Voltage localization in 
the superconductor case is illustrated in Fig. 2. In this 
case each bond has a random critical current, i^, and the 
whole network has the critical current, Ic- We solved the 
full fiow equations directly and then confirmed that at Ic 
voltage localizes on the minimum cut, i.e. the surface on 
which Ic = X^fceS *c is minimal. 

The easiest way to prove that random manifolds 
emerge at the macroscopic thresholds, Ic and Vc, is by 
using the cost function, v{i')di' . From this cost func- 
tion, it is evident that if cither the current or voltage 
can be kept to zero, then the cost itself is zero. In the 
case of a varistor, as the external voltage is increased, 
the network distributes the voltage drops so as to keep 
the current at zero. It succeeds in doing this until it is 
no longer possible and the first time this occurs is on the 
shortest path, occuring at the critical voltage, Vc- This 
argument applies to any cost function which has strictly 
zero current up to threshold. This means that the behav- 
ior of the local v — i characteristic after threshold has no 
influence on the value of Vc- In a similar manner, in the 
superconductor case, a large network attempts to keep 
the voltage at zero in all of the bonds until it is impos- 
sible to do so. As the applied current is increased, the 
current is distributed so as to keep the currents in all of 
the bonds below their threshold value, again in a very co- 
operative manner. The first applied current, Ic, at which 
voltage appears is determined by when a surface of bonds 
all have reached their local critical currents. This surface 
is the minimum cut. 

To make contact with experiment, note that varistors 
are materials which are insulating below a critical elec- 
tric field, Ec, after which they become highly conducting 
p^-p^ . The physics of their operation is understood to 
originate in the grain boundaries of polycrystalline ce- 
ramics. The interior of the grains of typical varistors 
are conducting, but the grain boundaries are insulating 
at low voltage (below about 3V). The onset of current 
flow occurs when the applied voltage is large enough to 
induce a path of grain boundaries to become conduct- 
ing. This physics is naturally encapsulated in lattice 
models in which a node represents a grain and a bond 
represents a grain boundary p^-p^] . However the grains 
have varying sizes and the grain boundaries have vary- 
ing properties, with some boundaries having higher onset 
voltages than others. That is, we must study disordered 
networks. These networks are also highly non-linear as 
has been demonstrated by studies of flow onset across in- 
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dividual grain boundaries ]lq ]. The linear dimension, L 
of the lattice models is related to the size of real varistors 
by L = l/g, where I is the sample size and g the grain 
size. Considerable progress has been made using numeri- 
cal simulations which have provided good agreement with 
experiment Here we make the connection with 

the statistical physics of random paths. 

As discussed above, network models of polycrystalline 
varistors have a current onset given by the minimum over 
all paths of Vc = J^kuP ■ This is just the shortest-path 
problem of Eq. (1), where the voltage threshold on a 
bond replaces the cost of that bond. The shortest-path 
problem is usually in the universality class of the directed 
polymer in a random medium(DPRM) |0 (execept for 
strong disorder ) so that Vc behaves like the energy 

of a DPRM, i.e. 

^4 = ali + a2i^ (3) 

where ai and 02 are independent of sample size and the 
exponent 9 is universal. Its value is known to be exactly 
1/3 for paths through two dimensional systems and 
0.248 ± 0.004 for paths through three dimensional me- 
dia 1^^. A key feature of (3) is that the threshold elec- 
tric field of the lattice Ec — Vc/ L is size independent, in 
contrast to the size effects produced by rare-fluctuation 
theories of electrical and dielectric failure |^J2^;|^ , where 
Ec approaches zero logarithmically in the large lattice 
limit. However the key difference is that in dielectric 
breakdown, local regions irreversibly make the transition 
to the conducting state, so that no voltage is required to 
maintain these regions in the conducting state after fail- 
ure. In contrast, high quality varistors are reversible so 
that there is a steady state current at a fixed applied volt- 
age. The path on which current flows is usually self-afhne 
and has length, Lp, given by P,^, Lp — biL + b^zL'', 
where 61 and 62 are dependent on the disorder distribu- 
tion, but the exponent C is universal. Its value is known 
to be exactly 2/3 [|^ for paths through two dimensional 
systems and 0.62 ± 0.01 for paths through three di- 
mensional media. The exponents 9 and C. are related by, 
9 = 2C — 1, which is consistent with the numerical re- 
sults quoted above. The difference in voltage between 
the lowest threshold path and the onset voltage of the 
next filamentary current path scales in the same way as 
the energy gap in the DPRM problem, which decreases 
logarithmically with increasing sample size |^^. In the 
strong disorder limit, the paths become highly tortuous 
and are no longer in the DPRM class p3| , p4| , p6[ . 

In the case of ceramic superconductors the low-angle 
grain boundaries have much higher critical currents than 
the high-angle grain boundaries ||2^. As in the varistor 
case, lattice models take the grain centers to be nodes and 
the links between grains to represent the grain boundaries 
1^,^. Grain boundaries may have a flux flow charac- 
ter or a Josephson Junction character or be resistive. If 
a grain boundary is resistive, its critical current is zero. 
As demonstrated above, the onset of voltage occurs on 



the minimum cut through such a network and this min- 
imum cut is related to domain walls in Ising magnets 
[p|,P HTl]j30[ | . The critical current of these networks then 
behaves in the same way as the energy of Ising domain 
walls, i.e., 

Ic - + 021" (4) 

where d is the spatial dimension, ci and C2 are indepen- 
dent of sample size, and = 1/3 in two dimensions and 
9 = 0.82 ±0.02 in three dimensions. There is thus a sam- 
ple size independent critical current density in these net- 
works, in contrast to the size effect which occurs in fuse 
networks Again this difference is due to the fact 

that these networks must maintain a current of at least 
ic on the A:"' bond in order for a voltage to appear there. 
The surface on which voltage localizes at Ic (i.e. the min- 
imum cut) is in the universality class of domain walls in 
the random-bond Ising model and so it is self-afhne and 
has asymptotic roughness, w, given by, w = C3L'', where 
C = 2/3 in two dimensions and C ^ 0.41 ± 0.01 in three 
dimensions |H,|ll],|o|. The exponents 9 and C are related 
to each other by the scaling relation, 6' = 2^ -I- d — 3. 

We have demonstrated that current or voltage localize 
on special manifolds in non-linear networks. As partic- 
ular examples, we showed that flow in a varistor begins 
on a path which is equivalent to a shortest path through 
the random medium, and in superconductors, voltage on- 
set occurs on a minimum cut through the network. The 
macroscopic current or voltage thresholds are then equiv- 
alent to the energy of a random manifold. Since these 
macroscopic thresholds only depend on the local current 
or voltage thresholds (and not on the local v ~ i charac- 
teristics after threshold) the manifolds on which current 
or voltage localize in varistors and superconductors are 
described by the universal directed polymer or random 
surface exponents. 
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FIG. 1. Behavior of a 50x50 varistor network, with the 
bonds having onset voltages uniformly distributed on the in- 
terval [0,2], asymptotic resistance one ohm, and with the ex- 
ternal voltage applied in the horizontal direction. The behav- 
ior of the average bond is the abrupt curve(smaller circles) 
in the top figure, while the behavior of the network is the 
smoother curve in the top figure(bigger circles). The localiza- 
tion of current on the shortest path is illustrated in the lower 
figure. This is the current pattern at V/ L = 0.59. At higher 
voltages further filamentary current paths emerge. 




FIG. 2. Behavior of a 50x50 Josephson junction network, 
with the bonds having onset currents uniformly distributed 
on the interval [0,2] and asymptotic resistance one ohm, and 
with the external current injected in the horizontal direction. 
The behavior of the average bond is the abrupt curve(smaller 
circles) in the top figure, while the behavior of the network 
is the smoother curve in the top figure(bigger circles). The 
localization of voltage on the minimum cut is illustrated in 
the lower figure. This is the voltage pattern at I/L — 0.58. 
At higher currents further sheet-like voltage surfaces emerge. 



